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Abstract 


In  this  thesis,  the  crystallographic  difference  Fourier 
method  is  employed  to  investigate  the  interaction  between 
the  enzyme  Glycogen  Phosphor ylase  a  and  some  of  its 
inhibitors.  Ihe  focus  is  on  the  "Nucleoside  site".  The 
interesting  finding  of  an  unusual  stacking  interaction 
between  the  enzyme  and  its  inhibitors  sheds  light  on  the 
specificity  of  this  site.  The  bi oenerge tics  of  the 
interaction  are  discussed  as  well  as  the  mechanism  of 
inhibition . 
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I.  Introduction. 


Glycogen  phosphor ylase  plays  a  major  role  in  various 
regulatory  processes  in  glycogen  metabolism.  It  degrades 
storage  glycogen  using  inorganic  phosphate  to  form 
glucose- 1 - phosphate  for  producing  AT?  through  the  glycolytic 
pathway.  The  enzyme  responds  to  energy  control  via  AM?  and 
AT?,  feed-back  inhibition  by  glucose  and  glucose-6 — 
phosphate,  and  a  number  of  hormonal  and  nervous  controls.  In 
addition  to  its  substrates  glycogen,  G-1-P  and  inorganic 
phosphate,  it  is  also  known  to  bind  nucleosides  and 
nucleotide  derivatives,  regulatory  enzymes,  divalent  metals 
and  glucose. 

A  substantial  body  of  research  has  been  done, 
especially  concerning  the  structural  and  functional  aspects 
of  the  enzyme.  The  enzyme  is  known  to  exist  in  two 
chemically  and  functionally  distinct  forms,  Phosphor yl ase  a 
(with  phosphate  on  serine-14)  and  Phosphorylase  b  (without 
phosphate) ,  the  interconversion  of  the  two  forms  is 
regulated  by  a  special  kinase  and  phosphatase. 

Phosphorylase  a  in  the  absence  of  ligands  is  partially  in 
the  active  (E)  form  whereas  phosphory lase  b  exists  mainly  in 
the  inactive  (T)  form.  The  B-I  equilibrium  is  regulated  by 
various  unknown  control  mechanisms,  using  either  metabolites 
or  macro  mo  1  ecu les  cr  both. 

Eoth  a  and  b  forms  of  the  enzyme  from  rabbit  muscle 
have  been  investigated  by  cry st alio graphic  studies  and  the 
structures  have  been  resolved  at  2.5  angstroms  and  3.0 
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angstroms  respectively  (Sprang  and  Fletterick  1979;  Weber  et 
al  1978) .  Distinct  binding  sites  lave  been  found  for  several 
of  the  inhibitors,  activators  and  substrates.  The  functions 
of  each  of  these  sites  and  their  effect  on  the  R-T 
conversion  are  listed  as  follows: 

Site/ligands  State  Functions 

AMP/ATP  R/T  Metabolic  response 

to  ATP/AMP  level 

Glucose/G- 1-P  T/R  Inactivation/  activation  of 

phosphorylase  a 
Nucleoside  1  Inactivation  of 

phosphorylase  a 

Glycogen  E*  Partial  activation,  dissoc¬ 

iation  of  tetramer  to  dimer 
Ser-r14/(p)  1/(5)  Hormonal  response 

To  understand  the  molecular  mechanism  of  R-T  conversion 
for  any  of  the  above  entries  (sites)  is  very  difficult.  One 
way  to  increase  our  understanding  is  to  lock  into  each  of 
these  sites  crystall cgra p bically  to  investigate  the 
interaction  between  effectors  and  the  binding  sites. 

The  purpose  of  this  thesis  is  to  study  the  nucleoside 
site  by  examining  a  number  of  fused  ring  compounds,  with 
respect  to  its  detailed  interactions  with  the  enzyme 
Phosphorylase  a. 

The  method  employed  is  known  as  difference  Fourier 
synthesis.  It  was  first  used  in  protein  crystallography  by 
Stryer  in  1964  to  examine  the  binding  of  azide  ion  to 


Myoglobin.  Since  then,  it  has  teen  extensively  used  ty 
protein  crystallcg rapher s  to  investigate  protein-ligand 
interactions,  protein  conformational  changes,  location  of 
heavy  atom  binding  sites  etc- 

The  difference  Fourier,  Delta  D,  is  defined  by 


a  D  =  y  I  m  |aF  1 L  *  &xj>  -  2.  l  A  X 

where  V  =  volume  of  unit  cell 

7 T\  =  figure  of  merit  (  a  measure  of 
the  correctness  of  the  phases, 

=1  represents  correct  phases). 

=  F ^  -  Fp 

Fp‘  =  structure  amplitude  of 
the  Protein  complex 

Fp  =  structure  amplitude  of  the  Parent 

=  phases  of  the  Parent  structure  amplitude 
The  power  of  the  method  lies  in  the  fact  that  a  large  number 
of  similar  structures  can  be  investigated  ty  using  only  the 
parent’s  phases  plus  the  individual  structure  amplitudes  for 
the  particular  binding  experiment.  The  resulting  density  map 
from  this  calculation  represents  the  difference  between  the 
protein  complex  and  the  parent.  In  the  case  of  ligand 
binding  the  peaks  represent  the  ligand  and  the  structural 
changes  induced  by  the  ligand  with  peak  heights  equal  to 
half  of  their  true  values.  A  detailed  decryption  of  the 


method  can  be  found  in  the  book  by  Blundell  and 


u 


Due  to  the  extensive  application  of  this  method,  it  is 
important  tc  look  into  its  reliability  and  accuracy.  This 
has  been  dealt  with  in  great  detail  by  Henderson  and 
Moffat  (1971).  In  their  study  they  focused  cn  errors  in  the 
difference  Fourier  maps  from  three  sources:- 

1)  use  of  the  difference  Fourier  coefficients  Delta  f 
instead  of  the  true  structure  amplitude; 

2)  experimental  errors  in  the  Delta  F  values; 

3)  experimental  errors  in  the  phases  of  the  parent 
structure  amplitude. 

The  remarkable  conclusion  they  have  drawn  is  that  the 
difference  Fourier  map  would  have  a  much  lower  error  level 
than  the  corresponding  Fourier  map  of  the  derivative 
structure  (which  would  reguire  a  new  complete  MIF.  phase 
calculation) . 
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A.  Materials  and  Methods- 

Cry st alii ration 

Single  crystals  of  Phosphor ylase  a  suitable  for 
diffraction  studies  vere  kindly  supplied  by  Dr.  Madsen  and 
S.Shechosky.  They  were  prepared  as  the  procedure  outlined 
below . 

Habbit  muscle  P hosphor ylase  b  was  prepared  by  the 
method  of  Fischer  and  Krebs  (1962)  with  manor  modifications 
including  the  use  of  5  mM-dithiothreitol  in  place  of 
mecaptoethanol .  Phosphorylase  a  was  prepared  from 
Phosphorylase  b  with  kinase  by  the  methods  of  Krebs  et  al 
(1964).  After  three  crystallizations  in  a  micro-needle  form 
from  10  raM  Sodium  glycerophosphate,  5  mM  Dithiothreitol, 

1  mM  EDT A  (pE  6.7),  the  Phosphorylase  a  was  transferred  into 
10  mM  3FS  (N , N-bis  (2-hydroxymethyl) -2  aminomethane  sulphonic 
acid),  0.1  mM  EDTA  (pH  6.7),  by  passage  through  a  Sephadex 
G25  column  at  39  degrees.  Portions  of  concentrated  stock 
solutions  of  magnesium  acetate,  dithiothreitol  and  glucose 
were  added  to  the  protein  solution  to  give  final 
concentrations  of  10,  3,  and  50  mM,  respec ti vely . 

The  protein,  at  concentrations  of  1.5  to  3  %,  was 
placed  in  small  glass  tubes  5  cm  long,  2.4  mm  internal 
diameter,  and  sealed  at  both  ends  with  Parafilm.  Although 
Phosphorylase  a  would  form  tetragonal  crystals 
spontaneously ,  a  small  portion  of  a  highly  diluted  stock 
suspension  of  finely  ground  seed  crystals  were  usually 
injected  to  control  the  number  of  nucleation  centers. 
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Crystallization  was  allowed  to  commence  at  22  degrees,  with 
storage  at  16  degrees  after  a  few  days.  This  crystallization 
procedure  has  been  reported  earlier. (Fletterick  et  al  1976). 

The  enzyme  crystallizes  in  tetragonal  space  group 
P43212  having  lattice  parameters  a=b=128.4  angstroms  and 
c=116.6  angstroms  with  one  monomer  of  97,400  daltons  in  the 
asymmetric  unit. 

Materials 

Caffeine,  FMN  Sodium  salt  and  Inosine  were  obtained 
from  Sigma  Chemical  Cc.,  ATP  Sodium  salt  was  obtained  from 
Terochem  Co.,  ATP  had  been  monitored  by  a  Paper 
Chromatographic  method  to  ensure  no  contamination  of  AMP. 
Ascending  chromatography  has  teen  used  as  described  by  the 
Biochemicals  reference  guide  (  PL  Biochemicals  Inc.  1977) . 
Euffer  system  No.  1  (  Isofcutjric  acid/conc. NH20H/  Water  in 

ratio  of  66/1/33)  was  used. 

Crystal  soaking 

The  ligand-binding  studies  that  are  referred  to  in  this 
report  were  carried  out  by  soaking  the  protein  crystals  in  a 
solution  of  standard  buffer  with  added  ligands.  The  standard 
buffer  is  10  mM  N, N- bis  (2-hydroxymethyl) 2-a minome thane 
sulphonic  acid  /  0.1  mM  FLT  A  /  10  mM  magnesium  acetate, 
pH  6.7.  The  soaking  time  and  ligand  concentration  are  listed 
in  Table  1. 

Data  Measurement 

All  data  measurements  were  done  on  a  modified  Syntex 
Analytical  instruments,  P21  automatic  diffractometer,  as 


J 


7 


TAELE  1 


Crystals  soaking  and 

other 

conditions 

(see  text) . 

Caff  eine 

Inosine 

FMN 

ATP 

No.  of  Crystals  used 

51 

9 

15 

14 

Total  no.  of  observations 

44009 

11165 

16809 

16398 

Total  no.  of  reflections 

21599 

10351 

11571 

126  15 

Total  overlaps 

37590 

1509 

8234 

6141 

Merging  R  factor 

0.  025 

0.035 

0.049 

0.022 

Delta  F/F 

14% 

5% 

7% 

6  % 

Radiation  decay  rate 

3.3% 

1-3% 

2.  1% 

2.9% 

(  percent  per  hour  ) 
Lattice  constants 

128.0 

128.2 

128.2 

128.2 

(in  Angstroms  ) 

116.4 

1  16.4 

116.5 

117.0 

Concentration  of  ligand 

2  m« 

1 5mM 

0. 5m M 

2  OmM 

reported  earlier  (Fletterick  et  al  1976)-  The  counter  and 
X-ray  source  were  positioned  55.0  cm  and  28.5  cm, 
respectively  from  the  crystal.  The  data  were  collected  using 
Ni  filtered  Cu  radiation  with  the  X-ray  tube  powered  at 
20001*.  The  diffractometer  was  operated  in  the  Omega-scan 
mode  with  7  steps  of  0.012  degrees  scanned  at  the  peak 
center.  From  these  7  measurements,  the  highest  neighboring  4 
were  selected  and  summed. 

The  background  correction  was  applied  by  analyzing  the 
background  variation  as  a  function  of  PHI  and  TWO  THETA.  CHI 
dependence  was  not  taken  into  account  since  the  beam 
cross-sect icn  is  circular.  Two  background  data  files  were 
collected  from  each  crystal,  the  first  consists  of  a 
measurement  of  background  intensity  as  a  function  of  PHI  at 
constant  ISO  THETA  (12.0  degrees)  and  the  ether  one  consists 
of  a  measurement  of  background  intensity  as  a  function  of 
TWO  THETA.  The  measurements  were  made  in  between  the  lattice 
points  by  doubling  the  lattice  parameters  and  measuring  a 
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selected  series  of  odd  values  of  h,k,l's.  Approximately  0.5 
degree  increments  of  TWO  THETA  and  10  degree  increments  of 
PHI  at  TWO  THETA  =  12.0  degrees  were  made,  andthe  Omega-scan 
technique  was  used  as  for  the  data  collection. 

Data  Processing 

Radiation  damage  was  monitored  by  repeatedly  measuring 
10  strong  reflections  uniformly  distributed  between  TWO 
THETA  of  3  tc  29  degrees,  and  analysed  as  a  function  of  time 
and  TWO  THETA  using  the  decay  model  of  Fietterick  et 
al(1976)  ,  Variation  as  function  of  h,k,l  was  not  been  taken 
into  account. 

Standard  Lorentz  and  polarization  corrections  were 
applied  as  well  as  absorption  correction  (North  and 
Phillips  1968)  which  was  analysed  as  a  function  of  PHI  only. 

Approximately  half  of  the  observable  HKLs  to  2.5 
angstrom  resolution  were  measured  for  the  difference  Fourier 
studies.  These  reflections  were  taken  from  lists  of 
reflections  that  had  been  sorted  according  to  the  product  of 
figure  of  merit,  magnitude  of  parent  structure  amplitude  and 
resolution  (Sin  THETA)  squared.  This  procedure  should  give 
the  best  measured  and  best  phased  data  with  a  constant 
distribution  with  respect  to  resolution.  Three  lists  were 
prepared  according  to  this  scheme  representing  reflections 
in  the  spherical  shells  of  4.5  -  20.0  angstroms,  3  - 
4.5  angstroms  and  2.5-3  angstroms  respectively.  The  top 
halves  of  these  three  lists  (strongest  F's)  were  used,  and 
from  them  HK1  data  lists  were  prepared  for  each  crystals. 
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About  15  crystals  were  used  to  complete  ore  data  set  because 
a  single  crystal  diffracts  for  about  14*r16  hours  before  it 
decays  severely.  Thirty  reflections  were  included  in  each 
data  list  for  scaling  together  the  data  sets  from  each 
crystal.  These  reflections  were  evenly  distributed 
throughout  the  reciprocal  lattice.  In  many  cases,  a  large 
degree  of  overlap  of  HKL • s  were  allowed  in  two  consecutive 
data  lists.  The  data  lists  were  merged  and  scaled  using 
G.N.Eeeke's  scaling  program.  The  merging  B  factor  for 
crystals  (Ec)  is  defined  as 

\ 


!  i 

SUM  SUM  |I  - I  | 
hkl  i  j  hkl  hkl  | 


Ec  = 


SUM  I 
hkl  hkl 

Where  Ihkl  is  the  intensity  after  a  weighted  scaling  of 
the  ith  measurement  of  the  intensity  of  reflection  hkl  and 
I (hkl)  is  the  mean  value  of  the  N  estimates  of  this 
intensity  from  different  crystals.  The  summation  is  taken 
over  all  equivalent  reflections  measured  for  each  crystals. 
Those  data  lists  which  showed  a  poor  E  value  (i.e.>10%)  were 
recollected.  The  final  E  values  for  each  data  set  are  listed 
in  Table  1. 

In  this  study,  5  different  data  sets  were  used,  namely 
Parent  (50mM  glucose)  ,  ATP,  Inosine,  FMN  and  Caffeine.  All 
of  these  were  at  2.5  angstroms  resolution.  Most  of  the  data 
for  the  Parent  and  Caffeine  were  collected  by  B.Fletterick 
and  ethers  (Sprang  and  Fletterick  1979). 
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Heavy  atom  derivatives,  ethylmercurithiosali- 
cylate  (EM1S)  ,  dimercury  acetate(DMA)  and  a 
"double- deri vative"  Pb-EMTS,  have  been  used  to  extend  the 
x-ray  diffraction  data  set  from  3.0  to  2.5  angstroms 
resolution  in  the  Multiple  Isomorphous  Replacement  (MIR) 
phasing  procedure  (Sprang  and  Fletterick  1979).  In  the  DMA 
experiments  caffeine  was  used  as  a  stabilizing  agent  since 
crystals  soaked  in  DMA  alone  showed  an  increased  radiation 
sensitivity  (from  6  to  7  %  loss  in  diffracted  intensity  per 
hour  at  3  angstroms  resolution) .  Inclusion  of  caffeine  at 
10  mM  in  the  soaking  buffer  reduced  the  radiation  damage  to 
3%  per  hour.  Caffeine  showed  up  as  the  fouth  highest  peak  in 
the  difference  electron  density  map  and  was  used  in  the  MIR 
phasing  procedure. 

The  difference  Fourier  synthesis  for  each  ligand 
binding  experiment  was  calculated  using  the  Delta  F  between 
protein  complex  and  parent  plus  the  most  recent  2.5  angstrom 
native  MIR  phases.  The  programs  used  are  DIFFER  and  FO DRIER. 
The  latter  performs  Fourier  synthesis  while  the  former 
calculates  Delta  F  and  scales  data  anisctropicall y.  This  is 
important  for  some  poorly  integrating  data,  but  it  is  not 
important  in  this  study  because  all  the  ligands  promote  the 
T  state  of  the  enzyme  which  is  already  present  in  this 
crystal  form,  so  the  data  are  isomorphic  with  respect  to  ail 
the  three  axes  X,  Y  and  Z. 

Graphics  System  (M M S-  X) 

The  MM5-X  (Molecular  Modelling  System)  graphics  system. 


purchased  from  Washington  University,  St. Louis,  was  designed 
primarily  for  crystallographic  model  building,  hut  at  the 
time  of  this  writing  is  a  research  tool  under  continual 
development  and  not  a  routine  working  instrument.  It  has  a 
large  number  of  built-in  options.  The  system  allows  the 
three  dimensional  presentation  and  interactive  manipulation 
of  molecular  structures  (i.e.  global  orientation  and 
translation  of  structures,  and  also  those  conformational 
changes  that  can  be  represented  by  bond  rotations  with  up  to 
six  bonds.),  together  with  background  displays  of  contoured 
electron  density  maps.  This  instrument  is  meant  to  be  a 
sophisticated  version  of  the  Richards*  box  discribed  below, 
del  fitting 

Preliminary  model  fitting  was  done  by  contouring  the 
different  density  maps  and  displaying  them  in  the  Richards* 
box  {Richards  1968).  Pig  1  represents  a  schematic  drawing  of 
the  Richard  box.  Since  the  scale  of  the  model  is  identical 
with  that  of  the  map,  the  corresponding  positions  in  the  map 
and  in  the  model  will  appear  superimposed  when  the  model  is 
fitted  correctly,  and  they  will  not  move  with  respect  to 
each  other  regardless  of  the  position  of  viewing.  Labguip 
Models  ( 1  cm/angstroa)  were  used  and  map  sections  were 
calculated  fcr  two  orientations,  one  perpendicular  tc  the 
crystallographic  z-axis  direction  and  the  other 
perpendicular  to  the  crystallographic  x-axis  direction. 

Cnee  the  rough  conformation  of  molecule  has  been  found, 
it  was  used  as  a  guide  for  choosing  the  most  appropriate 
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Fig.  1.  Schematic  drawing  of  the  Eichards  box. (see 

text) 

structure  from  published  papers  (of  small  molecules 
structures  which  best  represented  the  ligand;  Caffeine: 

Sutor  1  953;  Inosine:  Thewalt  st  al  1970  ;  F M N :  Eear  et 
al  1973;  A1P:  Kennard  et  al  19 '0),  which  is  then  fitted  into 
the  density  map  with  the  fIMS-X  system.  Co-ordinates  of  the 
ligands  were  displayed  on  the  screen.  By  doing  single  bond 
rotations,  global  rotation  and  translation,  the  models  were 
fitted  into  the  density  map  (which  was  contoured  in 
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3-dimensions  and  appeared  as  a  net  form)  to  give  the  present 
co-ordinates.  The  contour  level  for  these  maps  were 
arbitarily  assigned  to  provide  maximum  interpre tahilit y 
before  noise  level  gets  too  high  as  judged  by  number  of 
peaks  in  solvent  region. 

Stereo-drawings  of  the  fitted  ligands  were  done  by 
decoding  the  graphic  display  file  generated  by  the  MMS-X 
system  and  plottings  were  done  by  the  University  of  Alberta 
Calcomp  plotter. 

Co-ordinates  of  the  native  proteins  were  obtained  from 
the  model  building  in  this  laboratory  (Sprang  and  Fletterick 
1S79).  These  are  unrefined  and  contain  error  of  >  0.5 
Angstroms. 

Least  sguares  planes  through  the  atoms  of  the  bases 
were  calculated  by  the  method  of  Blow(1S60). 


Computer  Epograms 

The  authors  of  the  computer  programs  employed  in  the 
present  study  and  their  functions  are  listed  as  followed. 


NAME 

CEDNCHC 


CEY3 

DIFFER 


AUTHOR 

Jurgen  Svgusch 
George  Eeeke 
Eavid  Blow/others 


FCURIEF. 
PL  2 


George  Reeke 
Robert  Fletterick 


FUNCTION 

Data  reduction 

Data  merging 

Delta  F  calculation, 

scaling 

Fourier  synthesis 
Contour  electron 


density  map 


MAPPREP 

MMS-X 

STEREO 

ACCESS 

GRADIENT 


Stephen  Sprang 


Washington  U. 
Colin  Broughton 


T.  J.  Richmond 


Sprang  &  Fletterick 


Transfer  electron  density 
map  from  AMBAL  to  MMS-X 
Graphics  system 
Stereo  plot  of  MMS-X 
graphic  displays 
Calculate  accessible  sur¬ 
face  area  according  to  al¬ 
gorithm  of  Lee  &  Richards 
Calculate  electron  density 
gradient  at  given 


co-ordinates 


15 


B.  Jesuits. 

Koaenclatare  System 

She  rcmenclature  systems  for  the  molecules  studied  are 
presented  in  Jig  2.  Dihedral  angles  are  named  according  to 
the  system  cf  S undar alingam  as  illustrated  also  in  Fig  2. 


C  (U';  fc'Xo 


C(a')  Ea/DO 

C  C3'>  EXO 


Jig-  2.  Nomenclature  system  used  for  this  thesis.  The 
ribose  is  named  according  to  its  conformational  state.  If 
the  C2*  is  deviated  from  the  plane  (made  by  atoms  C4 * ,  C1s 
and  01* )  towards  the  same  direction  as  the  C5*,  the  ribose 
will  be  called  C2*  endo.  When  the  deviations  of  C2*  and  C3* 
from  the  plane  are  the  same,  the  ribose  can  be  called  either 
endo  or  exo,  whereas  when  the  deviation  is  dorminated  on  one 
side,  only  one  of  the  two  name  can  be  used  which  points  to 
the  atom  with  greater  deviation. 


J 
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Fig.  2 


Continued 
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Conformation  of  Molecules 

a  preliminary  fitting  was  carried  out  because  of  the 
ease  of  obtaining  the  overall  conformation  by  manual 
measurement  compared  to  the  MMS-X  system,  especially  when 
dealing  with  a  ribose  ring  as  in  Inosine  and  ATP  etc.  As 
mentioned  in  the  Materials  and  Methods,  in  the  preliminary 
model  fitting,  the  maps  were  contoured  twice  in  two 
orthogonal  directions  because  it  was  impossible  to  fit  the 
model  unambiguously  by  looking  from  one  direction  only. 
Contour  maps  in  the  direction  of  X-axis  placed  underneath 
the  model  provides  an  additional  guideline  for  the  fitting. 
Then  the  appropriate  structures  were  chosen  from  published 
papers  as  the  initial  structures  for  MMS-X  which  generated 
the  present  coordinates. 

Fig  3  (a-g)  illustrate  the  best  fit  of  the  Caffeine, 
Inosine,  FMN  and  ATP  molecules  to  the  difference  electron 


density  map  as  obtained  by  the  procedure  described  above. 
The  following  is  a  description  of  the  conformation  of  each 
individual  molecule. 


Caffeine : - 

Since  the  map  does  not  resolve  individual  atoms,  the 
caffeine  molecule  is  treated  as  a  rigid  body  when  fitted 
into  the  map.  Contouring  at  a  higher  level  of  electron 
density  should  indicate  the  position  of  the  ring. 
Considering  this,  we  have  obtained  two  possible 
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Pig  3(a-g)  .  These  figures  represent  the  interpretation  of 
the  difference  electron  density  map,  the  contour  level  for 
the  map  is  arbitarily  assigned  as  described  in  the  Materials 
and  Methods. 


Pig.  3a.  Caffeine  1 


Pig.  3b.  Caffeine  2 


1  9 


Fig.  3c.  Inosine 


Fig.  33.  FMN 
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Fig.  3e.  AT?  (SIN) 


Fig-  3f-  ATP 


(ANTI) 
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Fig.  3g.  ATP  in  the  SYN  conformation  with  a  higher 
contour  level,  of  electron  density. 

orientations  of  the  molecule  which  correspond  to  a  180 
degree  rotation  along  the  length  of  the  molecule.  It  is 
difficult  to  conclude  at  this  point  whether  one  of  the 
two  or  both  was  the  preferable  oreintation  based  on  the 
present  density  map.  Further  discussion  will  be 
presented  in  latter  section. 

I no sine: - 

During  the  preliminary  studies,  it  was  discovered  that 
the  inosine  molecule  exists  in  the  SYN-conf ormation 
which  is  different  from  the  preferred  ANTI-conf ormation 
in  solution  as  published  (Schweyer  et  al  1968;  Danyluk 
and  Hruska  1968) .  In  this  study,  only  the  C{2’)  endo  and 
C(3’)  endo-puckering  of  the  ribose  have  been  considered, 
since  these  are  known  to  be  the  most  stable  forms 
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(Sundaralingam  1969).  From  observation  in  the  Richards 
box,  it  is  found  that  only  the  C(2I)  endo-p ucker in g  of 
the  ribcse  ring  gives  a  reasonable  fit. 

The  dihedral  angle  CHI  of  the  molecule  is  232 
degrees,  which  falls  within  the  boundary  proposed  by 
Sundaralingam.  This  will  be  discussed  in  more  detail  in 
the  next  section. 

The  distance  between  the  0(5*)  of  the  ribose  and 
the  N (3)  of  the  base  is  found  to  be  2.6  angstroms. 
Therefore  the  formation  of  hydrogen  bond  between  these 
two  atoms  which  stabilizes  the  SYN-conf orma tion  is 
highly  likely. 


FMN:- 

Different  conf or maticns  of  the  rings  have  been  observed 
for  the  FMN  molecule  (Eear  et  al  1973;  Smith  et  al 
1977) .  Since  no  conclusion  could  be  drawn  on  the 
planarity  of  the  molecule  at  this  resolution  level,  the 
molecule  was  taken  as  copianar  in  this  study. 

Considering  the  fact  that  oxygen  is  more  electron 
dense  than  the  methyl  group,  the  part  of  the  map  bearing 
higher  electron  density  is  assigned  to  the  oxygen-- 
bearing  ring.  In  building  the  ribityl  tail  of  the 
molecule,  it  was  possible  to  arrange  all  the  atoms  in  a 
staggered  position  and  have  them  fitted  reasonably  well 
into  the  electron  density  map  at  the  same  time.  The 
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isoallcxazine  ring  can  be  fitted  vith  a  moderate  degree 
of  flexibility.  However,  due  to  the  presence  of  the 
ribityl  tall  this  flexibility  is  restricted.  The  density 
map  also  allows  a  limited  degree  of  freedom  for  the 
ribityl  tail,  but  these  variations  are  eliminated  if  the 
dihedral  angles  are  considered  to  be  restricted  in  a 
staggered  position.. 

The  portion  of  the  density  map  corresponding  to  the 
phosphate  group  is  net  well  defined.  This  may  be  due  to 
the  possibility  that  the  phosphate  group  which  lies  at 
the  solvent  region  of  the  protein  is  freely  rotating. 

Table  2  is  a  list  of  the  dihedral  angles  of  FMN , 

ATP  and  Jnosine. 

ATP:- 

Fitting  of  the  AT?  molecule  into  the  difference  density 
map  of  2.5  resolution  posed  a  rather  difficult  problem 
due  to  the  large  degree  of  freedom  of  the  AT?  molecule 
which  includes  the  large  number  of  dihedral  angles  and 
the  different  puckering  of  the  ribose  ring.  Since  the 
C(2')  endo  and  the  C{3')  endo  are  the  most  stable 
conformation,  only  these  two  types  of  sugar  conformers 
have  teen  considered  in  the  fitting  procedure.  When  the 
difference  density  map  was  contoured  at  the  level 
regularly  used  for  the  study  of  other  ligands,  it  was 
observed  that  the  contour  map  could  not  accomodate  the 
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Table  2 


Dihedral 

angles 

of  FMN , 

J.  no 

sine 

and  ATP. 

ATP 

X 

f 

9s 

uJ/ 

cj/ 

uJz 

246 

-  174 

130 

59 

160 

145  -32 

114 

Irosine 


E  "IN 


x  f 

232  ^75 


"XjL 

•*4 

X 

25 

176 

67 

160 

63 

164 

96 

whole  molecule.  Apparently  part  of  the  molecule  was  not 
resolved  at  this  contour  level.  Therefore,  another 
contour  tap  was  produced  at  a  much  lower  level  to 
accomodate  the  whole  molecule.  Taking  into  consideration 
that  although  the  phosphates  should  be  more  electron 
dense  than  the  rest  of  the  molecule,  they  are  likely  to 
be  moving  around  in  the  solvent  region,  then  the 
resulting  density  would  be  lower  than  that  of  the  ribose 
ring  plus  the  base.  Iherefore  the  denser  region  was 
assigned  to  the  ribcse  ring  and  the  base.  Both  the  C  (2 1 ) 
enao  and  C(3’)  endo  have  been  fitted  in  the  SYN  and  ANTI 
conformation.  It  was  clear  that  the  C(2f)  endo-SYN 
conformation  can  be  best  fitted.  However,  the 
possibility  of  having  an  ANTI  conformer  cannot  be 
excluded.  Since  AT?  strongly  favors  the  ANTI 
conformation  (Eerthod  and  Pullman  1973),  a  fitting  of 
the  map  with  an  ANTI  ATP  is  also  included  for  the  reader 
to  compare  the  two  fittings.  The  phosphates  can  be 
fitted  ir  a  number  of  positions  related  by  rotations 


. 


along  different  diphosphate  bonds  as  shown  in  fig3e.  As 
the  contouring  was  dene  at  an  intermediate  level#  only 


one  of  the  five  conformations  was  observable  which  was 
therefore  taken  as  the  most  preferable  position 
(fig  3  (e ,  f ,  g)  . 


In  studying  the  relative  conformations  of  the  rases 
of  these  ligands  with  respect  to  the  protein,  it  was 
observed  that  there  were  two  amino  acid  residues  in  the 
neighbourhood  of  the  ligands  which  sandwich  the  base  in 
a  stacking  fashion.  They  are  the  Tyrosine  612  and  the 
Phenylalanine  285.  fig  4  (a-e)  illustrate  the  stacking 
positions  for  Caffeine#  Inosine,  FUN  and  AT?-  They  are 
all  found  to  he  bound  to  the  protein  in  a  relatively 
coplanar  manner  except  FMN,  which  shows  a  deviation  of 
about  25  degrees  which  will  be  accounted  for  in  the 
later  section.  Moreover  all  these  ligands  adopt  a  more 
or  less  parallel  longitudinal  axis  when  stacked  (in 
ether  words,  if  the  direction  along  the  length  of  each 
molecule  is  defined  as  its  longitudinal  axis#  the  axes 
of  different  ligands  adopt  a  similar  orientation). 

Table  3  gives  the  equation  of  planes  of  the  various 
ligands.  The  angles  between  these  planes  have  also  been 
computed  and  listed  in  Table  4.  Co-ordinates  for  ligands 
are  listed  in  Table  5  (a-e)  . 

-rotein  -  Ligand  Interactions 

To  investigate  the  overall  structural  changes  of  the 
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Fig.  4(a-e)  are  illustrations  of  how  the  ligands  are 
between  the  two  amino  acids  ?he-285  and  lyr-612,  all 
ligands  are  being  projected  on  the  plane  parallel  to 
ring  of  Tyrosine  612. 


Fig.  4a.  Caffeine  1 


stacked 

the 

the 


Fig.  4b.  Caffeine  2 
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lig.  4c-  Inosine 


Fig-  Ud-  FMH 
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Pig.  4e.  AT P 


TAELS  3 


Equations  of 

least  square  planes 

of  the  ligands 

and  the 

two 

stacking  amino  acids,  XCcsA+YCos £+30050=0/  when 

e  D  = 

perpendicular 

distance 

from  the 

origin  of  the  unit  cell 

to 

the  plane  where  all  the 

atoms  of 

the  ring  lie. 

A ,  B  and 

C 

are  angles  between  the 

normal  to 

the  plane  and 

positive 

X, 

Y,  Z  axes. 

NAME 

CosA 

CosE 

CosC 

D 

?he-28  5 

-0.619 

0.655 

-0.433 

-16.53 

Tyr-6  1  2 

-0. 70S 

0.623 

-0. 329 

-10.98 

Caffeine  1 

-0.675 

0.687 

-0.269 

-  3.  88 

Caf f eine2 

-0.693 

0.653 

-0.307 

-11.  38 

Inosine 

-0.794 

0.545 

-0. 270 

-18. 72 

FMN 

-0.892 

0.313 

-0.328 

-32. 45 
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TAE1E  4 


Dihedral  angles  between  the 


NAME 

2he 

Tjr 

Car  f  1 

AT? 

12.6 

7.6 

2.  6 

EMN 

25.9 

20.7 

25.2 

Ino 

15.  1 

7.4 

10.6 

Caff  2 

cr 

• 

CO 

2.3 

3.1 

Caff  1 

10.  1 

5.4 

0.  C 

lyr 

8.8 

0-0 

Phe 

0.0 

2ast  squared  planes  of  table  2, 


Caf  f2 

XH2 

EMN 

AT? 

5.4 

1  1  .4 

26.2 

0.0 

22.8 

14.  9 

0.0 

8.7 

0.0 

0.0 

Table  5 (a-e)  are  lists  of  Cartesian  coordinates  of  the 
fitted  ligands  defined  in  the  unit  cell  with  the  standard 
origin  as  used  for  the  protein.  The  three  columns 
corresponds  to  the  x,  y  and  z  coordinates  in  Angstroms. 

TABLE  5a.  Caffeinel 


X 

Y 

Z 

N  1 

38.77 

36.64 

29.40 

C  2 

37.85 

35. 56 

28.92 

C  10 

39.70 

37.14 

28.36 

C6 

38.  79 

37.  18 

30.65 

N3 

37.03 

35.  15 

29.84 

Oil 

37.92 

35.  26 

27.  77 

C  12 

36.  06 

34.06 

29.46 

C  4 

36.  99 

35.67 

31.17 

C5 

37.  84 

36.64 

31.57 

NS 

36.27 

35.  22 

32.  1  7 

0  13 

39.64 

38.10 

30.90 

N  7 

37.  53 

36.  90 

32.92 

C  8 

36.60 

35.99 

33.22 

C  1  4 

38.  19 

37.  86 

33.  82 

Table  5b. 

Caf f eine2 

X 

Y 

Z 

N  1 

37.96 

35.75 

29.06 

C2 

38.  91 

36.89 

29.34 

CIO 

38.  05 

35.22 

27.69 

C6 

37.06 

35.21 

29.95 

N  3 

38.82 

37.32 

30.56 

Oil 

39.66 

37.  18 

28.  46 

C  12 

39.  72 

38.46 

30.93 

C4 

37.  93 

36.  81 

31.54 

C  5 

37.08 

35.80 

31.25 

N  9 

37.72 

37.28 

32.  74 

0  13 

36.34 

34.25 

29.55 

N  7 

36.  36 

35..  54 

32.43 

ce 

36.  77 

36.  50 

33.23 

C  14 

35.29 

34.54 

32.  64 

’  y ' 
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Tatle  5c.  Inosine 


X 


N  1 

37.  SO 

C2 

38.  71 

C6 

37.23 

N3 

38,  94 

C  4 

38.35 

C5 

37.  51 

06 

36.56 

N  7 

37.  06 

N9 

38.  45 

C8 

37.  65 

C  1 » 

39.28 

01  * 

38.  87 

C  2 f 

40.  70 

C  4  * 

40.  05 

02* 

41.28 

C  3  ' 

41.  15 

03* 

4  1.  20 

C5  * 

40.20 

05' 

40.68 

Y 

Z 

35.38 

29.32 

36.47 

29.20 

35.01 

30.51 

37.36 

30.12 

37.03 

31.31 

35.  94 

31.54 

34.01 

30.56 

35.96 

32.  86 

37.68 

32.  49 

37.  02 

33.37 

38.  88 

32.74 

39.77 

31.  67 

38.60 

32.53 

40.52 

31.26 

38.02 

33.  65 

40.02 

32.  19 

40.82 

33.33 

40.41 

29.81 

39.  12 

29.46 
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Table 


X 


N0 1 

38,  19 

N  03 

37.  IS 

CO  2 

37.97 

CO  2 

37.  75 

C04 

37.  00 

004 

36.45 

C04A 

37.45 

N  05 

37.31 

C  1 0  A 

38.  03 

N  10 

38.50 

C  05A 

37.65 

C06 

37.  43 

C  09  A 

38.  3  1 

COS 

38.80 

C  13 

39.27 

C  0  7 

37.  83 

C08 

38.  52 

C  1 2 

33.97 

C11 

37.59 

C  14 

39.  77 

0  14 

40.  7  1 

C  1  5 

40.59 

CIS 

40.  84 

C  16 

41.85 

0  16 

41.  59 

C  17 

42.37 

C  1 7 

42.76 

P  17 

43.  70 

018 

42.85 

0  19 

44.  77 

C20 

44.40 

.  FMN 

y 

Z 

37.  65 

31.03 

36.62 

32.94 

38.93 

32.  72 

37.81 

32.27 

35.47 

32.27 

34.55 

32.87 

35.  29 

30.88 

34.20 

30.31 

36.47 

30.31 

36.  35 

29.09 

34.07 

29.05 

32.80 

28.40 

35-  16 

28.39 

34.  98 

27.09 

37.55 

28.40 

32.70 

27.04 

33.71 

26.44 

33.55 

24.97 

31.35 

26.3  4 

38.27 

29.37 

37.51 

30.20 

39.38 

28.69 

40.38 

29.73 

3  8.96 

27.92 

37.75 

27.26 

39.96 

26.89 

41.  18 

27.  68 

42.23 

26.88 

43.31 

26.30 

42.  59 

27.85 

4  1.39 

25.74 
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Table  5e.  AT?  (SYN) 


X 


C2 

38.  72 

N  1 

37.86 

N3 

38.80 

C  4 

37.  94 

C5 

36.  97 

N  9 

37.81 

C  € 

36.93 

N  10 

36.  12 

N7 

36.  30 

C  8 

36.79 

C  1  1 

38.78 

C  2  * 

38.  75 

0  1 » 

39.96 

C  4  1 

41..  06 

02* 

38. OC 

C3  * 

40.  32 

03* 

40.47 

C  5  * 

41.67 

05  * 

42.  73 

?1 

42.  75 

0  1 

4  1.36 

011 

44.00 

06  * 

42.86 

22 

43.38 

02 

44. 92 

022 

43.  13 

07  » 

42.35 

?3 

40.  92 

03 

41.27 

033 

40.25 

0333 

40.32 

Y 

Z 

37.86 

33.00 

37.26 

33.91 

37.57 

31-73 

36.58 

31.38 

35.  96 

32.17 

36.03 

30.12 

36.33 

33.55 

35.81 

34.43 

35.00 

31.45 

35.07 

30.22 

36.  16 

2  8,.  9  6 

37.63 

28.51 

35.84 

29.59 

36.71 

28.97 

37.52 

27.27 

37.48 

28.  1  1 

36.86 

26.72 

37.  48 

30.  17 

38.55 

29.67 

39.88 

30.04 

40.46 

29.61 

40.  70 

29.  53 

3  9.86 

31.47 

40.94 

32.66 

40. 48 

32.92 

42.  27 

32.06 

4  0.55 

33.95 

40.  04 

34.  1  8 

38.53 

34.81 

40.83 

35.  1  5 

3  9.84 

32.90 
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protein  induced  by  these  ligands,  the  gradient  search 
program  was  used  as  described  in  Materials  and  Methods.  This 
program  calculates  the  electron  density  gradient  at  each  of 
the  atomic  coordinates  of  the  protein,  the  gradient  will  be 
greater  than  zero  when  there  are  structural  changes. 
Therefore  it  can  be  used  to  locate  xne  major  scrucrural 
changes  within  a  protein  if  there  are  any.  All  the  ligands 
have  been  examined  except  caffeine,  due  to  the  high  noise 
level  of  the  density  map  as  mentioned  previously.  No 
significant  main  chain  movement  was  observed,  which  agrees 
with  the  results  of  visual  inspection  of  the  density  maps 
and  the  small  value  obtained  from  the  %  structure  amplitude 
changes  as  listed  in  Table  1. 

In  checking  for  the  side  chain  movements,  only  the  area 
within  six  angstroms  from  the  surface  of  the  ligand  was 
studied.  No  significant  density  mass  was  observed  in  the 
neighbourhood  of  the  ligands,  so  it  can  be  interpreted  as 
absence  of  side  chain  movement  in  binding  these  ligands. 

The  inhibitor  site  contains  twelve  non-polar  uncharged 
amino  acids  (Pro281,  Phe2S5,  Phe2S6,  Leu37S,  Pro380,  Ala382, 
Leu383,  Pro610,  Met614,  Ala615,  Ile760  and  Met763),  Six 
polar  uncharged  amino  acids  (Asn133,  Asn282,  Asn28h,  Tyr572, 
Gly611  and  Tyr612),  six  basic  amino  acids  (Arg385,  His570, 
Lys573,  His613,  His767  and  Arg769)  and  six  acidic  amino 
acids  ( Asp283 ,  Glu2S7,  Glu381,  Glu38^,  Glu571  and  Asp759)  in 
a  box  within  4  angstroms  surrounding  the  ligands  (The  box  is 
defined  by  the  minimum  and  maximum  x,  y,  z  coordinates  of 


the  ligands  plus  an  extension  of  four  angstroms  along  each 
of  the  three  axes  in  both  directions.  It  appears  as  a 
rectangular  fcox)  .  Most  of  these  amino  acids  are  labelled  in 
Fig  5e. 

Table6a4b  are  lists  of  the  distances  between  the 
ligands  and  protein  side  chain  atoms  of  interest  (those 
atoms  pairs  which  may  form  hydrogen  bond  or  electrostatic 
interaction) .  When  the  box  dimension  is  extended  to  ten 
angstroms,  the  number  of  negative  charged  amino  acids  become 
slightly  larger  than  that  of  positive  charged  amino  acids.  A 
number  of  hydrogen  bonds  may  exist  between  AT?  and  the 
protein,  however,  due  to  the  uncertainty  of  the  atomic 
coordinates  in  this  model  plus  the  comparatively  lower 
quality  of  the  fitting  of  ATP,  a  firm  conclusion  has  not 
been  drawn  at  this  point.  Some  posible  pairs  are  N10  &  ND  of 
Asn282 ,  03*  4  022  of  Glu381  and  02’  4  0E1  of  Glu571.  It  is 
interesting  that  there  are  no  short  contacts  (less  than  four 
angstroms)  between  the  protein  and  the  phosphates  of  AT?. 
However,  the  are  four  glutamic  acids  and  one  arginine  in  a 
box  of  10  angstroms  around  the  05'  of  inosine.  Table6a  shows 
that  besides  the  stacking  amino  acids  ?he285  and  Tyr612,  the 
ligand  only  comes  into  contact  with  Asn282,  glu287,  His570, 
Glu57  1  and  Gly6  1  1 . 

As  mentioned  before,  all  the  ligands  are  sandwiched 
between  Tyr-612  and  ?he-285.  An  attempt  has  been  made  to 
investigate  if  these  two  amino  acids  have  been  tilted  due  to 
the  inclined  stacking  of  the  isoalloxa zine  ring  of  FMN.  No 
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TABLE  6a. 

List  of  all  van  der  Naals  contacts  between  the  ligands  and 

Protein  within  four  angstroms. 


CAFFEINE 

1 

DISTANCE 

HIS 

570 

C  1 2 

CE  1 

4.0 

GLY 

6  1  1 

N7 

0 

3.5 

C8 

0 

3.9 

C  1  4 

C 

3.8 

C  1  4 

0 

2.9 

CAFFEINE  2 


ASN 

282 

C  1  4 

ND 

3.5 

GLY 

611 

N9 

0 

3.4 

C8 

C 

3.8 

C3 

0 

3.4 

INOSINE 

GLY 

61  1 

N9 

0 

3.  7 

C  8 

c 

3.9 

C8 

0 

3.2 

Cl  1 

FMN 

0 

3.8 

ASN 

282 

004 

ND 

3.  3 

GLY 

611 

N03 

0 

3.6 

002 

c 

3.  6 

002 

0 

2.8 

CO  2 

0 

3.4 

ATP 


ASN 

282 

N  1 0 

ND 

3.2 

GLU 

287 

N  1 

OE2 

3.9 

N  1  0 

OE2 

3.8 

GLU 

38  1 

03' 

OE2 

3.5 

GLY 

611 

C2 

0 

3.6 

N  1 

C 

3.9 

N  1 

0 

3.0 

C6 

C 

4.0 

C6 

0 

3.5 

N  10 

0 

3.8 

GLU 

571 

02’ 

OE  1 

3.2 

■ 
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TABLE  6b. 

List  of  all  the  interesting  interactions  between  the  protein 
and  the  ribcse  plus  phosphates  {  that  is  possible  hydrogen 
bonding  pairs  and  eletrostatic  interacting  pairs  within  ten 

angstroms) . 


INOSINE 

DISTAi 

ASN 

282 

01  ' 

OD 

9. 9 

Cl  ' 

ND 

8.5 

02' 

OD 

8.4 

02  ' 

ND 

7.7 

C3  ' 

ND 

9.8 

GLU 

287 

0  1  ' 

OE  1 

9.1 

Cl  ' 

0E2 

7.  1 

02 ' 

0S1 

6.9 

02' 

CE2 

5.  1 

03' 

OE  1 

8.8 

03' 

CE2 

7.2 

05* 

0E2 

8.8 

LYS 

289 

02* 

NZ 

9.3 

GLU 

381 

01  * 

0E2 

9-4 

05* 

OE  1 

8.4 

05 ' 

0E2 

6.  8 

GLU 

571 

01  ' 

OE  1 

7.4 

01  ' 

0E2 

9.4 

C3  * 

OE  1 

9.9 

05' 

CE 1 

6.2 

05  ' 

OS  1 

7.  9 

ASP 

768 

05' 

CD2 

9.9 

AEG 

769 

05* 

NE 

9.7 

05* 

NH2 

7.8 

05  ' 

NH  1 

8.1 

FMN 

ASN 

282 

004 

ND 

3.3 

014 

OD 

9.7 

014 

ND 

8.6 

GLU 

287 

0  14 

OE  1 

9.7 

014 

0E2 

7.6 

0  15 

OE  2 

9.3 

GLU 

38  1 

0  1  4 

CE  1 

8.0 

0  14 

OE  2 

6.7 

015 

OE  1 

9.4 

015 

0E2 

7.7 
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016 

QE  1 

5-  7 

C 1  6 

0E2 

4.  1 

017 

CE  1 

8,7 

017 

0E2 

6.9 

217 

0E  1 

9.  3 

P  1  7 

0E2 

7.6 

013 

0E2 

8.2 

019 

0E2 

S .  6 

020 

OE 1 

8.3 

020 

0E2 

6.7 

GLU 

571 

014 

OE  1 

7.0 

C  14 

0E2 

8.8 

015 

OE  1 

6.7 

0  15 

0S2 

8.3 

016 

OE  1 

C  C 

Jo  ^ 

016 

0S2 

6.8 

017 

CE  1 

7  a  0 

017 

0E2 

7.9 

P 1  7 

OE  1 

7.9 

PI  7 

0E2 

8.5 

C  1  8 

OE  1 

7.6 

0  18 

0E2 

8.  C 

019 

OE  1 

9.3 

019 

0E2 

9.9 

020 

CE  1 

8.0 

020 

0E2 

8.3 

AS? 

759 

?  1  7 

OD  1 

9.9 

018 

OD  1 

8.5 

018 

QD2 

8.  9 

ASP 

768 

016 

OD  1 

9.3 

016 

0D2 

7.  8 

017 

OD  1 

8.7 

C 17 

CD2 

8.0 

PI  7 

OD  1 

7.9 

P 1  7 

0D2 

7.  5 

018 

OD  1 

7.6 

013 

0D2 

7.4 

0  19 

OD  1 

9.  C 

C  1 9 

0D2 

8.  7 

020 

0  D  1 

6 . 8 

020 

0D2 

6.3 

AS? 

769 

C  1  8 

NE 

9.4 

C  1 8 

NH2 

8.6 

0  18 

NH  1 

7.4 

019 

NK2 

9.3 

019 

NH  1 

8.2 

020 

NE 

8.0 

020 

NH2 

7.1 

020 

NH  1 

5.  8 

i 
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ATP 

ASN 

282 

N 10 

ND 

3-2 

01  * 

OD 

9  .  C 

01  ' 

ND 

7.  9 

02' 

ND 

9.8 

?3 

OD 

9.8 

P3 

ND 

8.9 

03 

OD 

8-  6 

C3 

ND 

7.9 

033 

ND 

9.1 

0333 

OD 

8.8 

0333 

ND 

8.  7 

GLU 

287 

0  1  • 

OE  1 

9.8 

01  * 

0E2 

7.  7 

05  ’ 

0E2 

9.1 

PI 

0E2 

9.4 

0  1 

OE  2 

9.5 

06  ' 

0E2 

8.  4 

P2 

0E2 

8.6 

02 

CE2 

9.2 

022 

0E2 

S .  7 

07  ' 

0  E 1 

8.7 

07' 

0E2 

7.2 

?3 

OE  1 

7.7 

P3 

OE  2 

6.1 

03 

0  E  1 

6.  5 

03 

0E2 

4.9 

C33 

CE  1 

7.5 

033 

OE  2 

6.1 

0333 

0  E  1 

8.  3 

0333 

0E2 

6.5 

N  1 

0E2 

3.9 

N  1  0 

0E2 

3.8 

IYS 

289 

07' 

NZ 

9.6 

GLO 

381 

oi~ 

OE  1 

7.  1 

01 » 

0E2 

6 . 0 

02' 

CE  1 

7.3 

02  1 

0E2 

5.6 

03’ 

OE  1 

5.  2 

03’ 

0E2 

3.5 

05' 

OE  1 

7.9 

05  ' 

0S2 

6.6 

PI 

0  E 1 

9.0 

PI 

0E2 

7.6 

Cl 

CE  1 

9.3 

01 

0E2 

7.6 

01  1 

OE  1 

9.4 

Oil 

0E2 

7.9 

06’ 

CE2 

8.8 
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GLO 
01  * 
01  ' 
02* 
02* 
03  ' 
03’ 
05* 
05* 
PI 
PI 
01 
01 
011 
06’ 
033 
AS? 
02* 
03* 
03* 
05  * 
Cl 
Oil 
AEG 
01  * 
01  ’ 
Cl  • 
02  ’ 
02* 
02 » 
03' 
03  ' 
03’ 
05' 
05’ 
05' 
PI 
El 
01 
0  1 
01  1 
Oil 
06* 
06  ’ 


571 

OE  1 
0E2 
0  E  1 
CE2 
OE  1 
0E2 
OE  1 
CE2 
OE  1 
0E2 
OE  1 
0E2 
OE  1 
OE  1 
0  E  1 

768 
0D2 
OD  1 
0D2 
OD  2 
0D2 
0D2 

769 
NS 
NH2 
NH  1 
NE 
NH2 
NH  1 
NE 
NH  2 
NH  1 
NE 
NH2 
NH  1 
NH2 
NH  1 
NH2 
NH  1 
NH2 
NH  1 
NH2 
NH  1 


6*6 

8,  5 

3.2 

5.3 

4.7 

6.0 

7.8 
9-2 
8.1 

9.  5 

6.9 

8.4 

8. 9 

9.2 

9.1 

9.  2 

9.5 

7.8 

9.9 
9.9 
9.8 

9.5 

7.3 

8.4 

8.1 
6.2 

7.3 

6.6 

4.5 
5.  4 

9.7 

7.8 

7.9 
8.  7 

8.6 
8.6 
8.5 
9.0 

8 . 4 

9.9 
9.9 
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TAELE  7 

Surface  area  buried  and  the  energy  involved  (see  text) .  Due 
to  limitation  of  space  -abbreviations  are  used  for  this 
table.  When  the  name  of  the  ligand  first  appears  it  is  given 
the  full  name ,  XXX-C  =  ligands  -  protein  complex,  XXX  +  C  = 
ligands  plus  protein.  XXX-C  is  a  measured  quantity  while 
XXX+C  is  a  calculated  guantitj.  Delta  Gh  =  hydrophobic  free 
energy.  Delta  3tr  =  translational  and  rotational  free 
energy.  Delta  Gd  =  free  energy  of  dissociation  and  Delta  H  = 
enthalpy  of  the  interaction. 


Surface 

C  bange 

in 

T  A  S 

Name 

area 

area 

AGh 

AGtr 

Ajd  A  h 

Protein 

4536.8 

Caf f eine 1 

348.8 

Caff  1  +  C 

4885.6 

Caff  1-C 

4409.6 

476.0 

12.0 

-18 

-5.5  -11.5 

Caffeine  2 

351.3 

Caf f2+C 

4888. 1 

Caf f 2-C 

4409 . 7 

478.4 

12.0 

-18 

-5.5  -11.5 

FMN 

618.6 

FMN+C 

5155.4 

FMN-C 

4501.8 

653.6 

16.3 

-2  0 

-6.8  -10.5 

Inosine 

382.7 

Ino+C 

4919.5 

Ino-C 

4508.9 

410.6 

10.3 

-18.5 

-4.0  -12.2 

conclusion 

could  be 

drawn  due 

to  the 

ambiguity 

of  the  map. 

Sur Jace  area  calcula tion 

To  understand  and  interpret  the  ligand  -  protein 
interaction  is  very  difficult.  Recently  Janin  and  Chothia 
attempted  to  explain  this  problem  by  using  the  concept  of 
hv drophobicity  (Janin  and  Chothia  1978).  Changes  in 
accessible  surface  area  as  defined  by  Lee  and 
Richards  (1971)  had  been  used  to  determined  the  degree  of 
hy drophobicity  change.  Using  the  same  concept  the  surface 
area  changes  of  the  ligands  before  and  after  binding  to  the 
protein  have  been  calculated  and  listed  in  Table  7.  This 
will  be  described  below  in  more  detail. 


• 
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Energy  for  the  Binding  of  Ligands 

laking  the  empirical  value  of  25  cal  per  sguare 
angstrom  of  accessible  surface  area  buried  (Chothia  1S74), 
the  hydrophobic  free  energy  has  been  calculated  as  listed  in 
Table  7.  The  number  of  25  has  been  derived  from  experiments 
on  the  solubilities  of  amino  acids  in  different  solvents  by 
measuring  the  free  energy  of  transfer  of  amino  acid  from  a 
ncn-polar  reference  solvent  tc  an  aqueous  medium.  Fig  5&6 
are  space  filling  and  stick-bond  drawings  which  provide  a 
better  visualization  of  the  actual  ligand-protein 
interaction. 


C.  Discussion. 

3 a 5 - X  System 

Since  aas-X  has  been  used  throughout  the  study  for  the 
interpretation  of  the  results,  it  is  desirable  at  this  point 
to  comment  cn  seme  of  the  features  that  make  it  superior  to 
the  Richards  Eox  system.  Firstly,  since  both  the  molecule 
and  the  density  map  can  be  rotated  freely  and  can  be 
observed  from  any  orientation  desired,  it  gives  a  much 
better  fit  than  the  Richards  box.  Secondly,  since  the  system 
produces  the  measurement  of  coordinates  automatically,  this 
avoids  the  random  errors  generated  by  manual  measurement  as 
well  as  the  systematic  error  due  to  the  Richard  box  such  as 


' 
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Figures  5  & 
of  the  four 
ligands  are 


6 (a-d)  are  sticX-bond  and  space  filling  drawings 
ligands  being  studied,  the  orientation  of  the 
all  the  sane. 


Fig.  5a.  Caffeine 


Fig.  ca.  Caffeine 
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Fig.  6b.  Inosine 
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Fig.  5c.  FUN 


Fig.  €c.  FMN 
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Pig.  5d.  ATP  (SIN) 


fig .  6d .  ATP  (SIN) 
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lig-  5e.  ATP  (SIN) .  All  amino  acids  are  labelled  with 
one  letter  code. 

the  diffraction  (con traction)  of  the  half-way  mirror  and  the 
photo  distortion  due  to  xeroxing  of  the  contour  maps. 
Thirdly,  the  contour  level  of  the  map  can  he  changed  freely 
and  figures  can  he  prepared  more  easily  and  fitting  can  be 
criticized  more  readily. 

Stacking  interaction . 

The  stacking  of  the  ligands  between  the  Tyr-612  and 
?he-285,  is  a  complex  T)  T\  (pi-pi)  ,  van  der  tfaals 
interaction.  With  the  exception  of  some  electrostatic 
interactions  with  separation  in  the  order  of  nine  angstroms, 
no  hydrogen  bonds,  charge  or  salt  bridges  have  been 
observed,  though  hydrogen  bonding  might  be  suspected  between 
the  ring  nitrogen  of  Histidine-570  and  06  of  Inosine. 


-/  ■ 


However,  since  tie  distance  between  these  two  atoms  is  4.9 
angstroms,  the  formation  of  hydrogen  bonding  is  highly 
unlikely  even  allowing  for  the  errors  in  the  coordinates. 

Since  stacking  is  the  only  significant  force 
stabilizing  the  interaction  between  the  protein  and  ligands 
it  is  expected  that  the  degree  of  specificity  of  the  bindin 
is  limited  to  the  overall  molecular  structure  of  the  ligand 
This  observation  is  consistent  with  the  results  of  kinetic 
studies  of  rabbit  muscle  enzyme  which  show  that  the  enzyme 
has  relaxed  specificity  for  a  number  of  compounds  of 

ring  structures  ranging  from  purine  analogues  to 
iscallcxazine  ring  of  riboflavin  (Kasvinsky  et  al  1978b, c). 

The  present  study  is  not  the  first  incidence  of 
co-enzyme  stacked  between  two  protein  side  chains.  Similar 
binding  has  been  reported  of  FMN  in  Flavodoxin  (Watenpaugh 
et  al  1973),  in  which  case  the  stacking  involved  parallel 
stacking  with  a  tyrosine  and  45  degree  inclined  stacking 
with  a  tryptophan. 

Stacking  has  been  regarded  as  an  important  factor  in 
stabilizing  polynucleotide  structures  and  in  directing  the 
binding  of  aromatic  compounds  to  nucleic  acids.  In  a  review 
by  Bugg,  the  base  stacking  pattern  of  over  seventy  crystal 
structures  of  nucleic  acid  constituents  in  polynucleotides 
had  been  examined  (Bugg  et  al  1971).  It  was  found  that 
extensive  overlapping  between  bases  is  unusual. 

Surprisingly,  results  from  theoretical  calculation  also 
indicate  that  the  highest  geometric  overlap  between  bases 


f' 
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does  cot  cecessarily  lead  to  the  most  favorable  stacking 
interaction  (Orstein  et  a_l  1978).  In  the  present  study,  as 
indicated  by  fig  4a-e,  the  ligands  overlap  extensively  with 
Phe-285  but  show  only  partial  overlapping  with  Tyr-612.  In 
the  case  of  AI2,  the  adenine  base  does  not  penetrate  deeply 
into  the  binding  site,  so  the  base  ring  is  not  stacked  as 
extensively  as  the  other  ligands. 

In  an  attempt  to  locate  the  stacking  energy  minima 
between  bases,  Orstein  has  calculated  the  potential  energy 
of  stacking  as  a  function  of  the  rotational  angle  and 
separation  distance  by  decomposing  the  total  base^pair 
stacking  energy  into  its  components,  i.e.  electrostatic, 
polarization,  dispersion,  and  repulsion.  It  was  observed 
that  the  electrostatic  component  controls  the  location  of 
stacking  minima,  the  most  favorable  relative  nase-pair 
orientation.  Since  the  presence  of  heteroa toms  in  an 
aromatic  hydrocarbon  results  in  an  asymmetric  eleccr ostatic 
environment,  this  is  likely  to  control  the  mutual 
oreintaticn  of  the  bases.  As  the  number  of  heteroatoms  and 
their  location  in  the  ring  vary,  in  the  absence  of  specific 
hydrogen  bonding,  the  mutual  orientation  between  the  ligands 
and  the  stacking  amino  acid  side  chains  differ  with 
different  ligands.  In  addition,  the  remaining  parts  of  the 
molecule  also  impose  a  restriction  on  the  orientation  of  the 
base,  because  they  tend  to  position  themselves  in  a 


relatively  open  area  which  is  the  solvent  region. 

As  mentioned  in  the  Results,  the  ligands  were  stacked 


. 
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in  a 


more  or  less  parallel  fashion  with  respect  to  Iyr-612 
and  Phe-285,  with  the  exception  of  FUN,  which  exhibits  a 
plane  dihedral  angle  of  25  degrees.  In  the  stacking  cf  most 
bases,  a  plane  dihedral  angle  of  0  ■“  5  degrees  is  guite 
common.  The  fact  that  the  FNN  has  such  a  large  inclinition 
may  be  accounted  for  by  its  increased  size  due  to  an 
additional  ring  which  causes  the  molecule  to  tilt  slightly 
in  order  to  accomodate  itself  in  the  binding  site.  If  it 
were  not  tilted,  it  would  collide  with  His-570. 

^22^0 tmaticn  of  nucleotides.. 

The  aly cosidic  torsional  angle  CHI  is  used  to  indicate 
the  relative  orientation  of  the  ribose  with  respect  to  the 
base  in  the  nucleotide.  It  is  defined  by  the  rotation  of  the 
C8-NS  bond  of  the  base  about  the  N9-C  1  *  bond  with  respect  to 
the  CI'-OI'  bond  of  the  sugar.  The  torsional  angle  is 
considered  to  be  zero  degrees  when  C8-N9  are  cis-planar  with 
respect  to  C1»-01«.  A  positive  CHI  value  is  obtained  as  the 
base  rotates  in  a  clockwise  direction  when  viewed  along  the 
C1-N9  axis.  The  ribose  is  in  the  ANTI  conformation  when  CHI 
is  equal  to  0+/-9Q  and  SYN  when  CHI  is  equal  to  180+/-90 
degrees  (Fig .  2) . 

In  the  present  study,  it  is  observed  that  both  AT  I  and 
inosine  are  most  likely  in  the  SIN  conformation,  with 
dihedral  angles  being  246  and  232,  respectively .  In  a  review 
by  Sundaralingam  investigating  crystal  structures  of  25 
nucleotides  and  nucleosides,  he  observed  three  distinct 
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regions  for  the  CHI  values  of  BETA  Purine  nucleoside, 
namely,  3-55,  120-  123  and  210-258  degrees.  In  the  last 
angular  range,  the  conformation  has  been  observed  to  be 
stabilized  by  05-H.-.N3  hydrogen  bonding,  and  all  the 
nboses  exhibit  C(2’)  en  dop  uckering  (  Sundaralingam 
1969,  1  975)  . 

Our  scudy  with  inosine,  which  exhibits  similar 
characterist ics  would  well  fit  into  the  last  group.  ATP, 
however,  due  to  the  absence  of  hydrogen  bonding  between  05 
and  N 3 ,  can  not  be  compar ed  wi th  these  gro ups  of  s tr uc t ures 
investigated.  From  another  study  by  Sundaralingam  it  was 
also  observed  that  nucleotides  are  conf ormaticnally  more 
’’rigid"  thar  nucleosides  (Butin  et  al  1972).  In  addition,  it 
was  proposed  that  nucleotides  exhibit  only  ANTI  conformation 
about  the  glycosyl  bond,  while  nucleosides  may  exist  in  both 
SYN  or  ANTI  conformation.  The  applicabili ty  of  the  arove 
generalization,  however,  would  be  limited  by  the  size  of  the 
sample  studied,  which  is  restricted  by  the  number  of 
structures  resolved  by  X-ray  crystallography.  Since  it  is 
questionable  whether  Sundaralingam * s  proposal,  which  was 
based  entirely  on  the  survey  of  X-ray  crystallographic 
results  (in  solid  state  in  crystals),  represents  a  genuine 
fundamental  difference  with  respect  to  the  properties  of  the 
two  groups  of  nucleo-compounds,  a  theoretical  investigation 
which  employs  the  PCILO  computations  has  been  performed  by 
Pullman  and  associates  (Eerthod  and  Pullman  1973a, b).  From 
the  study  of  a  number  of  purine  and  pyrimidine 
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5  nucleotides,  it  was  demonstrated  that  although  guaccsine 
ard  inosine  5  ’-phosphates  show  a  preference  for  the  SYN 
conformation,  most  nucleotides  prefer  the  ANTI  conformation. 

Pullman's  results  appear  to  lend  some  support  to  the 
previous  proposal  as  well  as  to  provide  a  better  theoretical 
consideration  of  the  problem.  However,  one  has  to  be  aware 
of  the  limitation  of  the  methodology  of  most  theoretical 
calculations.  In  this  case,  the  pseudo— rotation  of  the 
ribose  has  not  been  dealt  with  due  to  the  complexity  of  the 
problem.  The  ribose  was  prefixed  at  either  the  C  ( 2  * )  endo  or 
the  C  (3 1 )  endo  conformations.  Therefore  the  applicability  of 
the  results  of  the  calculations  is  questionable. 

Although  nucleotides  bound  to  enzymes  are  dikely 
distorted  from  their  preferred  conformations  in  solution, 
the  glycosyl  torsion  is  mcst  stubborn  to  distortion 
( Sundaralin qam  1975).  There  has  not  been  any  report  of  5* 
Adenine  nucleotide  with  the  SYN  conformatioi.  In  the  present 
study  the  ATP  was  observed  to  be  in  the  less  preferred  SYN 
conformation.  The  contradiction  between  this  observation 
with  the  cited  generalization  can  be  accounted  for  by  the 
following 

1)  The  above  generalization  may  not  be  representative, 

2)  There  might  be  some  factors  operating  in  this  particular 
case  which  may  give  sufficient  energy  to  overcome  the 
glycosyl  torsion  energy  barrier. 

The  first  possibility  has  been  considered  in  the  above 
discussion.  However,  it  is  interesting  to  note  than  ATP  when 


-r 


52 


bound  to  another  site  (AMF/ATP  activator  site)  of  the  same 
enzyme  also  exhibits  an  ANTI  conformation  (Yang  et  al 
unpublished  result  1378),  which  gives  another  incident  to 
support  the  validity  of  the  generalization  being  discussed. 
Therefore,  we  have  to  go  into  the  second  suggestion  to  find 
a  possible  explanation  of  this  unusual  binding  conformation. 
As  has  been  pointed  cut,  the  molecules  has  to  obtain 
sufficient  energy  a)  to  change  from  its  "preferred"  ANTI 
conformation  in  solution  to  the  SYN  conformation  in  the 
binding  site,  b)  to  have  enough  energy  to  stabilize  the 
"high-energy"  SYN  conformation  in  the  binding  site.  Ihe 
second  energy  requirement  can  be  easily  met  by  the  stacking 
energy  generated  after  the  ligand  has  been  bound  to  the 
protein.  The  question  now  to  be  considered  is  therefore  how 
the  molecule  undergoes  its  ANTI  -  SYN  transformation.  Let  us 
consider  a  hypothetical  situation  in  which  ATP  exists  in  the 
solvent  in  high  molarity,  there  is  finite  probability  that 
some  molecules  may  have  acquired  sufficient  energy  for  the 
ANTI  -  SYN  transformation.  Since  only  the  SYN  conformer  can 
be  fitted  into  the  density  map,  it  is  likely  that  that  the 
enzyme  can  stabilize  this  less  stable  conformation.  As  has 
been  pointed  out,  this  is  possible  only  in  the  case  of  a 
high  concentration  of  ligands. 

The  above  idea  has  been  corroborated  experimentally. 
Kinetic  studies  carried  cut  by  Kasvinsky  et  al  (1378b, c) 
demonstrated  that  with  adenine  and  hy poxanthine,  the 
addition  of  a  ribose  does  not  alter  the  inhibitory  constant 
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of  the  free  base  (Table  8).  However,  addition  cf  a  phosphate 
group  alters  the  inhibitory  constants  of  these  ligands, 
especially  in^  the  case  of  adenosine.  The  change  of  free 
energy  in  binding  from  incsine  to  lap  and  from  adenosine  to 
AMP  are  0.2  and  1.6  Kcal  respectively.  It  is  important  to 
look  into  the  reason  for  this  significant  difference  between 
the  energy  requirement  for  the  binding  of  the  two  molecules, 
.following  the  argument  cited  previously  in  this  section,  we 
can  assume  that  1)  although  in  solution  both  inosine  and 
adenosine  favor  an  ANTI  conformation,  these  nucleoside 
structures  are  rather  free  in  converting  from  one 
conformation  to  the  other,  2)  addition  of  a  phosphate  group, 
however,  fixes  the  nuclec-compcunds  in  their  preferred 
conf ormations  which  is  ANTI  for  AM?  and  possibly  SYN  for 
IMP,  3)  the  binding  site  of  these  ligands  favors  stacking  of 
SYN-conf ormers. 

It  can  then  be  postulated  that  in  order  for  AM?  to  be 
bound,  it  has  to  acquire  a  certain  amount  of  energy  to 
convert  to  a  rather  unstable  SYN  conformation,  whereas  in 
the  case  of  IMP,  this  is  unnecessary.  The  slight  increase  of 
the  inhibitory  constant  of  IMP  with  respect  to  inosine  may 
be  due  to  the  elec trosta tic  interaction  between  the 
phosphate  and  protein  side  chains  {possibly  glu287,  glu381, 
glu5 7 1  and  arg769,  tabled),  the  electrostatic  interaction 
may  be  expressed  by  the  eguat ion (Schulz  and  Schirmer  1879) 
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TABLE  8 

Kinetic  data  for  the  ligands  studied 
(Kasvinsky  et  al  1978b, c) . 


Liqands 

C  on  c  .  m  M  » 

Ki  mM. 

Delta  G b — & 

Caffeine 

1.0 

0.  1 

-5 1  5 

FMN 

1.0 

0.0  1 

-6.8 

Hypoxan thine 

0.9 

1. 1 

-4 . 1 

Inosine 

1.0 

1. 1 

-4 . 1 

IMP 

1.5 

-3.9 

Adenine 

1.0 

0.  4 

-4.7 

Adenosine 

1.0 

0.  4 

-4 . 7 

AMP 

— 

6.0 

-3.0 

ATP 

— 

60  (10-20) 

(-2.5) 

E  = 


3  3  2  % 

6.  R,x 


Where 


fc- 

e  = 
E  = 


partial  charges  in 
electron  charges  e  . 
distance  between  partial 
charges  in  Angstroms, 
dielectric  constant 
repulsion  or  attraction  energy 


in  Real  per  angstrom 


if  the  dielectric  constant  is  taken  to  be  30,  the  repulsion 
between  two  unit  negative  charges  at  distances  of  3,  6,  9, 
12  and  15  angstroms  will  be  3.7,  1.8,  1.2,  0.9  and  0.7  Real 
respectively.  From  table6  it  is  evident  that  there  are  four 
negatively  charged  amino  acids  (glu287,  glu381,  glu571  and 
asp768)  and  one  positively  charged  amino  acid  (arg769) 
within  10  angstroms  of  05'  of  inosine  in  the  oinding  site, 
if  a  phosphate  is  attached  to  the  05*,  the  net  force  can  be 
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as  high  as  4  times  1.2  Kcal  minus  0,9  Kcal,  that  is  3.S 
Kcal,  therefore  the  contrib ut ion  of  electrostatic  forces  to 
the  binding  of  IMP  would  be  significant.  This  value  of  3.9 
Kcal  for  a  mcnopole  approximation  may  account  for  the  fact 
that  AT?  does  not  go  into  the  binding  site  far  enough  and 
hence  resulted  in  a  weaker  binding.  However,  due  to  the 
large  number  of  charges  involved,  detailed  calculation  of 
the  net  force  is  difficult  and  has  not  teen  undertaken. 

In  the  case  of  ATP,  there  is  a  variation  in  the 
reported  Ki  values  between  solution  and  crystal  studies. 
Taking  the  results  of  crystal  studies  and  assuming  the  mean 
to  be  15  mM ,  we  can  show  that  the  free  energy  reguired  for 
the  binding  of  ATP  to  the  enzyme  is  2.5  Kcal  (Table  8). 
Therefore,  difference  in  free  energy  of  binding  between  AT? 
&  adenosine  is  2.1  Kcal  ,  this  may  be  the  results  of  1) 
electrcsta tic  interaction,  and/or  2)  conformational 
transition.  Further,  the  difference  between  binding  of  AT? 
and  AM P  indicates  that  the  presence  of  two  additinal 
phosphates  does  not  affect  the  binding  of  the  nucleotide  to 
the  enzyme  tc  any  great  extent  once  the  molecule  has  been 
made  ‘'rigid'1  by  the  first  phosphate. 

Bio-energetics  of  ligand-Protein  interaction 

Since  Kauzmann  introduced  the  concept  of  hydrophobic 
bonding  (Kauzmann  195S),  it  has  been  considered  by  some  to 
be  the  most  important  force  in  protein  folding,  formation  of 
protein  multisubunit  complexes. 


and  binding  of  small 
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ligands.  Using  the  concept  of  accessible  surface  area 
defined  by  lee  and  Richards  (1971)  as  illustrated  in  Fig. 7, 
it  was  proposed  that  energy  is  released  when  protein 
accessible  surface  is  “buried”,  that  is,  removed  from 
contact  with  solvent;  and  the  magnitude  of  this  energy  has 
been  determined  empirically  as  25  cal  per  sq.  angstrom  of 
protein  accessible  surface  area  buried  (Chothia  1974). 

Janin  and  Chothia  have  recently  applied  this  concept  to 
study  the  binding  of  coenzymes  to  Flavodoxin  (Janin  and 
Chothia  1979),  Lactate  dehydrogenase  and 
Glyceraldehyde-3-phosphate  dehydrogenase.  They  calculated 
Delta  G  from  the  following  eg uat ions : - 
Delta  G=  Delta  Grt  +  Delta  Gh 
Delta  G  =  RT  InKi 

Where  Delta  Grt  is  the  free  energy  loss  in  translational  and 
rotational  degrees  of  freedom.  Delta  Gh  is  the  gain  in 
entropy  due  to  hydrophobic  interaction  and  Ki  is  the 
inhibitory  constant. 

The  first  equation  is  derived  from  the  classical 
thermodynamic  equation  Delta  G  =  Delta  H  -  T  Delta  S.  The 
term  T  Delta  S  has  been  replaced  by  two  components  Delta  Grt 
and  Delta  Gh.  The  first  one  gives  a  negative  contrib ution 
while  the  latter  gives  a  positive  contribution  to  the 
entropy  of  the  system.  The  enthalpy  term  Delta  H  in  the 
equation  has  teen  assumed  to  he  negligible,  that  is,  no 
conformational  changes  or  other  energy  interactions.  For 
although  hydrogen  bonding  and  Van  der  Kaals  contact  do  exist 
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Pig.  7.  This  is  to  give  an  illustration  of  the 
accessible  surface  area  as  defined  by  lee  and  Sicharas  ( 1 9 7 1 ) 

between  the  protein  and  ligands,  they  would  make  similar 
interaction  with  the  solvent  when  not  bound,  making  that 
balance  negligible.  The  value  of  Delta  G  as  calculated  from 
the  two  equations  were  comparable,  hence  justifying  the 
omission  of  Delta  H. 

The  same  calculation  has  been  carried  cut  in  our  study 
using  FHN,  irosine  and  caffeine.  large  discrepancies  have 
been  found  in  the  Delta  G  calculated  from  both  equations 
(Table  7).  The  reason  for  the  seemingly  contradictar y 
finding  between  the  two  studies  can  be  due  to  the  following 
reasons.  1)  The  ligands  studied  by  us  were  not  buried  as 
completely  as  those  in  the  previous  study.  Hence  the 
hydrophobic  interaction  is  not  large  enough  to  compensate 
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for  the  entropy  less  of  the  co-enzymes  involved.  2)  Stacking 
plays  an  important  role  in  the  protein-ligand  interaction 
energy  in  our  study  while  rhis  does  not  exist  and  so  has  not 
been  considered  in  Janin’s  study.  This  stacking  interaction 
is  an  enthalpy  driven  process  and  can  not  be  accounted  for 
by  the  hydrophobic  energy  alone.  Therefore  an  appreciable 
Delta  H  term  must  be  considered  in  the  present  study. 

A  simple  rearrangement  of  the  equations  cited  above 
would  give  the  followings:- 
-Delta  H  =  ETln  Ki  -  T  Delta  S 

=  ETln  Ki  -  Delta  Grt  -  Delta  Gh 
If  this  Delta  H  is  contributed  solely  by  the  stacking 
interaction.  Delta  H/2  would  give  us  the  stacking  enthalpy, 
since  the  stacking  interaction  in  our  study  involved  2 
pair-stacking,  which  is  tyrosine  -  ligand  and  ligand  - 
phenylalanine . 

The  experimental  and  theoretical  value  of  Delta  H  for 
stacking  has  been  under  extensive  investigation  by  other 
workers  (Ts'o  et  al  1963;  Vorst  et  al  1965;  Stoesser  et  al 
1966;  Gill  et  al  1967;  Gupta  et  al  1378) .  No  agreement  has 
been  arrived  at  due  to  the  difficulties  involved  in  its 
measurement.  A  number  of  values  have  been  reported  ranging 
from  3  Kcal/mole  to  10  Kcal/mole.  The  Delta  H/2  values 
obtained  from  our  study  have  been  listed  in  Table  7.  They 
fall  well  within  the  reported  range  of  values  and  agree 
quite  closely  with  one  another.  This  suggests  that  the 
stacking  energy  involved  in  the  binding  of  these  ligands 
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does  not  differ  appreciably  from  ore  case  to  another.  The 
difference  of  their  binding  constant  is  therefore 
contributed  mainly  by  the  difference  in  their  entropy, 
which,  in  this  case,  is  the  rotational  translational  entropy 
and  the  hydrophobic  entropy.  Since  Delta  Gh  represents  the 
hydrophobic  interaction  between  the  ligand  and  the  protein, 
the  above  observation  suggests  that  the  difference  in  the  Ki 
of  one  system  from  another  is  determined  mainly  by  the 
difference  in  the  surface  area  of  the  ligand  buried  in  the 
protein . 

Mechanism  of  inhibition  of  enzymatic  f unction . 

From  the  above  discussion  it  has  gradually  become 
apparent  that  stacking  interaction  plays  an  important  role 
in  selectivity  of  the  nucleoside  binding  site.  As  mentioned 
in  the  introduction,  the  enzyme  phosphorylase  a  exists  in  a 
state  of  I-E  equilibrium  in  solution.  In  a  recent 
crystallographic  study  with  the  enzyme,  it  was  found  that  as 
the  enzyme  undergoes  its  1-S  conversion  induced  by  UEPG,  one 
of  its  polypeptide  loops  which  joins  helix  265-275  and 
290-310  is  displaced  from  its  "original”  position  as 
demonstrated  in  Fig  8.  A  close  examination  of  a  wire  model 
suggests  that  as  long  as  the  "movable"  loop  remains  in  its 
"original"  position,  the  entrance  of  substrate  into  the 
active  site  would  be  blocked  and  hence  the  enzyme  would  be 
in  its  inactive  form.  In  addition,  the  center  of  the 
pyrophosphate  of  UDPG  is  2.0  angstroms  from  the  GG  of 
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fig.  8.  Interpretation  of  the  difference  density  map 
from  UPPG  experiment,  red  represents  the  negative  contouring 
while  blue  represents  the  positive  contouring.  Usually  when 
the  chain  movement  is  will  defined,  one  can  correlate  the 
negative  density  with  the  positive  density,  however,  in  this 
case,  the  movable  loop  which  is  shown  in  dashes  might  likely 
becomes  disordered  so  no  positive  density  can  be  assigned  to 
it  as  its  final  position. 


Asp-283  of  the  (T)  form  which  would  suggest  that  the  Asp-283 
must  move  upon  binding  of  UDPG,  resulting  in  the  movement  of 
the  whole  loop. 

One  of  the  amino  acids  Phe-285  which  is  involved  in  the 
stacking  complex  is  located  in  the  movable  loop.  The 
presence  of  Inhibitors  attaches  this  Phe  to  Tyr-612  to  form 
a  stacking  complex,  thus  fixing  the  "movable"  loop  to  its 
"original"  position.  As  a  result,  the  I-form  of  the  enzyme 
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is  stabilized.  The  argument  is  supported  by  the  finding  that 
in  Potato  ph csphorylase ,  which  exhibits  no  regulation,  the 
?he-28  5  is  replaced  by  Serine  (Nakano  et  al  1978).  However, 
this  may  not  be  the  only  mechanism  contributing  to  this 
stability.  It  has  been  reported  that  this  inhibition  site  is 
Synergistic  with  respect  to  glucose  inhibition  of 
phosph-ory lase  a.  In  previous  study,  it  has  been  observed 
that  hydrogen  rending  may  exist  between  the  C2  of  the 
glucose  and  the  ND  of  Asp-284,  which  is  also  located  on  the 
movable  loop  (Yang  et  al  1978  unpublished  results) .  It  can 
be  postulated  that  this  hydrogen  bonding  may  give  further 
stability  to  the  loop  and  hence  the  enzyme  in  its  T  form.  If 
this  postulation  is  true,  this  can  provide  a  possible 
explanation  to  the  above  mentioned  Synergism. 

It  has  been  demonstrated  that  in  liver  phosphor ylase  a, 
a  functional  "nucleoside  sire"  also  exists,  both  in  vitro 
and  in  vivo.  This  "nucleoside  site"  appears  to  have  a  higher 
selectivity  than  the  muscle  enzyme,  for  it  only  allows  the 
binding  of  purine  and  purine  analogues.  However,  it  is 
questionable  whether  a  physiologica 1  effector  does  exist 
which  is  able  to  bind  to  this  "nucleoside  site"  in  vivo 
(Kasvinsky  et  al  19  78c)  .  To  answer  this  question,  a  number 
of  experiments  have  been  done  in  this  laboratory  by  soaking 
crystals  of  muscle  phosphcrylase  a  in  liver  extract  (Madsen 
and  Fletterick  unpublished).  Subsequent  crystallographic 
studies  have  revealed  an  unidentified  ligand  (s)  which  binds 
to  this  "nucleoside  site".  The  electron  density  map,  as 
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shown  in  Jig  9  suggests  that  the  hound  ligand  is  likely  to 
he  a  purine  analogue.  Since  the  muscle  enzyme  used  has  a 
relaxed  specificity  on  finding  of  effectors,  in  is  possible 
that  this  density  map  may  represent  the  average  of  a  number 
ligands  of  similar  structures.  Therefore  any  attempt  to 
interpret  this  density  map  has  to  be  done  with  caution. 
However,  since  the  site  may  be  involved  in  the  insulin 
governed  glucose  homeostasis,  it  is  expected  that  future 
identification  of  the  effector  in  the  physiological  system 

bring  about  significant  insight  on  the  understanding  of 
this  important  hormonal  regulatory  system. 


-ig.  5.  Unknown  ligands  in  the  liver  extract 
text)  . 
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D.  Conclusion 

In  the  present  study,  a  numner  of  fused-ring  ligands 
have  been  used  to  investigate  one  of  the  regulatory  sites  of 
the  muscle  enzyme  phosphorylas e  a.  Although  the  study  has 
been  somewhat  limited  by  the  resolution  of  the  difference 
electron  density  map  due  to  the  present  aifi  phases 
attainable,  we  were  able  to  obtain  a  number  of  interesting 
findings  on  the  energetics  of  the  protein-ligands 
interaction  as  well  as  some  aspects  on  the  inhibition 
mechanism.  However  we  have  not  been  able  to  obtain  a  full 
understanding  of  the  pattern  involved  in  the  stacking.  This 
may  he  rendered  possible  in  the  future  after  the  protein 
phases  have  been  refined  or  with  higher  resolution  Delta  ? 
maps  of  these  and  additional  ligands. 
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